The size-extensitivity of correlation energy estimators 
based on effective characteristic polynomials ^ 



Herbert H. H. HomeierS 



Institut fiir Physikalische und Theoretische Chemie, Universitdt Regensburg, 
D- 93040 Regensburg, Germany 



Estimators Hn for the correlation energy can be computed as roots of 
effective cliaracteristic polynomials of degree n. The coefficients of these 
polynomials are derived from the terms of the perturbation series of the 
energy. Prom a fourth-order M0ller-Plesset (MP4) calculation one can 
calculate with negligible effort a size-extensive estimator 112 that is in 
many cases much closer to the full CI correlation energy of the ground 
state than the MP4 value. [H.H.H. Homeier, J. Mol. Struct. (Theochem) 
366, 161 (1996)] Here, we prove that the estimators Hn for n > 2 are 
size-extensive if they are calculated from the MP series. 
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1 Introduction 



Many-body perturbation theory is a convenient tool to estimate the correla- 
tion energy of molecular systems. Usually, one calculates an estimate for the 
correlation energy by term-by-term summation of the M0ller-Plesset (MP) 
series 

E = Eo + Ei + E2 + E3 + Ei + Er, + ... , (1) 
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partial sum 

= x: (2) 



that is usually denoted by MPn. But this approximate value does not opti- 
mally exploit the information content of the terms Ej of the MP series. Better 
estimates can be obtained by using convergence acceleration or extrapolation 
methods to sum the perturbation series, as for instance Fade approximants or 
methods based on effective characteristic polynomials. The method of effec- 
tive characteristic polynomials has been introduced by Cizek and coworkers 
recently. It has been applied to the summation of the divergent perturba- 
tion series of anharmonic oscillators [0] and to correlation energies of model 
systems [Q-Q. Extensions of the method for the simultaneous treatment of 
several perturbation series have recently been proposed but this will not 
be considered in the following. 



If only perturbation energies up to the fourth order are available then one can 



use the 112 estimator |[T^. The 112 estimator is obtained as a root of a second- 



degree effective characteristic polynomial. The coefficients of this polynomial 
are related to the terms of the perturbation series. The 112 estimator can 
be calculated easily from the terms Eq, . . . ^E^ of the MP series and is size- 
extensive. As is well-known, the latter property is important for the treatment 
of larger systems. For a series of benchmark systems the 112 estimator proved 
to be relatively accurate as compared to a number of other estimators pO|] . 



Analogous estimators Iln can be derived from effective characteristic polyno- 
mials of higher degrees n. In the present contribution, we sketch the method 
of effective characteristic polynomials, basically to fix notation. Then, it is 
proved that the estimators Iln for n > 2 are also size-extensive if the under- 
lying perturbation theory is size-extensive, as is the case for the MP series. 



2 Correlation energy estimators based on effective characteristic 
polynomials 



In this section, we sketch the method of the effective characteristic polynomi- 
als. The characteristic polynomial Pn{E) of degree n in the unknown energy 
E has in the linear variation method the form 

Pn{E) = det \{<P,\H\<Pk)-ES,^k\ (3) 
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where are n orthonormal basis functions, and H is the Hamiltonian. 

If the Hamiltonian can be written as H = Hq + f3V, the polynomial has the 
form (i, Eq. (3.2)) 

n n-j 

Pn{E) = Y.E= E/".>/5' (4) 

j=0 k=0 



with fn,nfi = 1- This leaves L + 1 coefficients fn,j,k to be determined where 
L = n{n + 3)/2 — 1. In the method of the characteristic polynomial, they are 
obtained from the coefficients Ej of the perturbation series for E 

oo 

E = Y.E,P^. (5) 

j=0 

Since Pn{E) = for an eigenvalue E, one demands 

Pn{Eo + (3Ei + (3^E2 + ...) = 0(/3^+^) . (6) 



This means that the first L + 1 coeffients of the Taylor series (in /3) of the 
left-hand side of this equation, i.e. up to the coefficient of all have to van- 
ish. This yields a linear equation system for the unknown fn,j,k- We assume in 
the following that this linear equation system yields a unique solution for the 
coefficients fn,j,k- After the determination of the coefficients, the effective char- 
acteristic polynomial is fixed and also denoted by Pn{E) = Pn[Eo, . . . , El]{E) 
in order to make the dependence on the terms Ej explicit. The effective charac- 
teristic equation Pn[Eo, . . . , El]{E) = may then be solved for E. The lowest 
root is called Iln[Eo, . . . , El] or, more simply, Iln, if the values of the Ej are 
plain from the context. 

In the case n = 2 one obtains for /5 = 1 the simple expression 



E-}, — E-i 



El 

n2 = + El + — 

2 E2 -£'4 



Ei 



+ 



El 



{E2 — E^Y — 4 {E2 E4 



El) 



1/2 



{E2 -E4 — E. 



2\2 
3) 



(7) 



3 Proof of the size-extensitivity 



Consider a supersystem composed of M identical, non-interacting subsystems. 
Then, the true energy of the supersystem is the M-fold of the energy of a single 
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subsystem. An approximate method for the computation of the energy is called 
size-extensive if the approximate energy for the supersystem is the M-fold of 
the approximate energy of a single subsystem. 

It is well-known that MP perturbation theory is size-extensive order by order. 
This means that for all A; > the k-th order term of the perturbation series 
for the supersystem equals M Ek if the k-th order term of the perturbation 
series of a single subsystem is Ek- This implies that the MFk estimator is 
size-extensive for each order k. 

In order to prove size-extensitivity of the Ilr;, estimators with n > 2, one thus 
has to show that for each n and M the equation 

Un[M Eo, . . . , M El] = MUn[Eo, . . . , El] (8) 
holds where, as before, L = n{n + 3)/2 — 1. 

For given M and n and for given Eq, . . . ,El and /5, we can consider the fn,j,k 
and hence, also the effective characteristic polynomial 

n n—j 

P^[Eo, ...,El]{E) = Y^ E^ fn,j,kP'' (9) 

j=0 k=0 

as known. We introduce a new polynomial by 

P„(^) =M"P„[Eo, . . .,El]{E/M) (10) 

n n-j 

= E^' E[fn„kM--^]l3' ■ (11) 
j=0 k=0 

Eq. (|T0|) defines the polynomial, while Eq. { ^A\ ) shows that Pn{E) can also 
be regarded as an effective characteristic polynomial with new coefficients 
fn,j,k = fn,j,kM'^~K Note that fn,n,Q = 1 holds as required. 

Since 



Pn{M Eo + (3 M El + [3^M E2 + . . .) 

= M"P„[Eo, . . .,EL]{Eo + f3E, + (3^E2 + . . .) 

= , (12) 

the polynomial Pn{E) is identical to the effective characteristic poljTiomial for 
the energies M Eq, . . . , M El, i.e. 

M" P^[Eo, El]{E/M) = P^[MEo, ...,MEl]{E). (13) 



4 



Thus, the complete pattern of roots is scaled by M, since if e is any root 
of P„[^o, • • • , EL]iE) then Me is a root of P„[MEo, . . . , MEi](E). But this 
proves Eq. (H) since Iln[EQ, . . . , El] is a root of Pn[Eo, . . . , El]{E). 



4 Concluding Remarks 

Size-extensitivity is an important property that is helpful for obtaining reli- 
able energy estimates for larger systems. It is hoped that correlation energy 
estimates on the basis of effective characteristic polynomials become more 
widespread since they offer a conceptually and computationally relatively sim- 
ple but accurate computational tool as demonstrated in Only the MP 
series has to be computed, and its terms be combined to obtain an effective 
characteristic polynomial and the correlation energy as one of its roots as in 
Eq. ([^. Some error control is possible by comparison to some other methods 
for convergence acceleration of the perturbation series [ p!0[] . 

An extension of the method was proposed by the author in ||^. It allows to 
combine information from the terms of several short perturbation expansions 
for a small number of states. These simultaneous perturbation series can also 
be used for the construction of a characteristic polynomial. In this way, the 
somewhat demanding relation between the degree of the characteristic polyno- 
mial and the maximal order of the perturbation calculation is largely avoided. 
Applications of this extended method and a proof of its size-extensitivity are 
currently under investigation. 
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